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Lattice thermal conductivities of zincblende- and wurtzite-type compounds with 33 combina¬ 
tions of elements are calculated with the single-mode relaxation-time approximation and linearized 
phonon Boltzmann equation from first-principles anharmonic lattice dynamics calculations. In 9 
zincblende-type compounds, distributions of phonon linewidths (inverse phonon lifetimes) are dis¬ 
cussed in detail. The phonon linewidths vary non-smoothly with respect to wave vector, which is 
explained from the imaginary parts of the self energies. It is presented that detailed combination 
of phonon-phonon interaction strength and three phonon selection rule is critically important to 
determine phonon lifetime for these compounds. This indicates difficulty to predict phonon life¬ 
time quantitatively without anharmonic force constants. However it is shown that joint density of 
states weighted by phonon numbers, which is calculated only from harmonic force constants, can be 
potentially used for a screening of the lattice thermal conductivity of materials. 


I. INTRODUCTION 


At temperatures near or above room temperature, the 
thermal conductivity of non-metallic solids is governed 
by phonon-phonon collisions, which is an anharmonic 
phenomenon. Recent progress of computers, numerical 
algorithms, and methods of first-principles calculations 
have enabled us to predict anharmonic force constants 
quantitatively, and it has even become possible to repro¬ 
duce well the lattice thermal conductivities of compounds 
using those anharmonic force const ant s.EHUl However a 
detailed analysis of anharmonicity which is at work to 
produce the lattice thermal resistivity is still not well in¬ 
vestigated. 

In this work we performed systematic lattice thermal 
conductivity calculations for a class of compounds with 
the zincblende- and wurtzite-type structures for 33 differ¬ 
ent combinations of chemical elements. These two crys¬ 
tal structures are equivalent in local coordinations and 
different in their stacking orders. This is similar to the 
difference between the face-centered cubic and hexagonal 
close-packed structures. 

We mainly focused on the phonon lifetime since we 
believe that it is the most critical physical property to 
seek for a mechanism of lattice thermal conduction. In 
fact the lattice thermal conductivity and phonon life¬ 
times are related straightforwardly (see Sec. IIBI in 
the single-mode relaxation-time (SMRT) method, being 
merely proportional to each others. The results of these 
calculations are presented in Fig. [T] and Fig. [5] 

In Fig. Eit is found that the differences of the two crys¬ 
tal structures, i.e., the difference in their stacking orders, 
affects little the lattice thermal conductivities. The same 
behavior is also observed for the phonon linewidth distri¬ 
butions and phonon density of states which are presented 
in Appendix [G] For the ease of presentation this allows 
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FIG. 1. (color online) Lattice thermal conductivities calcu¬ 
lated at 300K with respect to volume per formula unit. 


us to focus our discussions only on the zincblende-type 
compounds. Moreover, since the zincblende-type com¬ 
pound has cubic symmetry, the phonon linewidth dis¬ 
tribution in the Brillouin zone will be considered to be 
nearly isotropic. Therefore the phonon linewidth distri¬ 
butions will mostly be investigated along T-L direction. 

Fig- [2] shows that a good agreement is obtained with 
experiments, and evidence the lifetime as an important, 
controlling, parameter. The lattice thermal conductiv¬ 
ities have also be computed from a full solution of the 
linearized phonon Boltzmann equation (LBTE)P As pre¬ 
sented in Appendix |T| the numerical results, with and 
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FIG. 2. (color online) Comparison in lattice thermal conduc¬ 
tivities at room temperature between experiments and calcu¬ 
lations with RTA and witout isotope scattering. 


without isotope scattering, are not very much different 
from those presented in Fig. [T] and Fig. [2] This con¬ 
firms our claim that the anharmonic lifetime should be 
understood carefully. Effect of structural disorder, which 
can vary phonon frequency, linewidth, and lattice ther¬ 
mal conductivity“tl^ was not considered in this study 
since the effect is expected to be small for most of these 
compounds. 

The paper is organized as follow. In section |II A| the 
method of calculation is reviewed, with the most techni¬ 
cal aspects described in Appendix [X]to[E| The results are 
discussed in section [TTT1 which is separated in two parts. 
In the first part we study the dynamical aspects of the 
scattering from the (w, q) maps of the self energy, and 
in the second part the average value of the lifetime over 
the Brillouin zone. The numerical results are collected in 
Appendix [F] to |T] and compared to experiments in Table 

hyi 


II. METHOD OF CALCULATION 

A. Phonon lifetime 


The phonon lifetime that we discuss in this study is 
computed from the imaginary part of phonon self en¬ 
ergy. The crystal potential is expanded with respect to 
atomic displacements. The coefficients of the third order 
which contain the anliarmonicity are used to compute the 
imaginary part of the self energy. For the Hamiltonian 


Jif the expansion can be written as 

J4? = 4> 0 + ST + 3^2 + ^3 + • • • (1) 

with 

& [u a (lK)] 2 , ( 2 ) 

Ikol 

^2 E K')u a (lK)up(l '*/), (3) 

Iko. Vk'{5 

=1J2 E E i"k") 

Ikol l' k'( 3 l"K ,,r y 

X U a (ln)'Ui3(l'K')u-y(l” k”), (4) 

where 4> 0 , and J^ n correspond to the constant poten¬ 
tial, kinetic energy, and n-body crystal potential terms, 
respectively, u (Ik) is the atomic displacement of the k- 
th atom in the l- th unit cell, m K is the atomic mass of 
type k. and a, /3, and 7 are the Cartesian indices. 
and denote the harmonic and cubic anharmonic 

force constants, respectively. The harmonic Hamiltonian 
is defined as =z£h = + ^2 ■ 

An atomic displacement operator may be written as^ 1 ^ 


u a (ln) = 
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( 5 ) 


where N is the number of unit cells in the crystal and h 
is the reduced Planck constant, r (Ik) is the equilibrium 
atomic position and and a^j are the phonon creation 
and annihilation operators of the normal mode of band 
index j and wave vector q. The harmonic frequency 
and polarization vector W (k, qj) are obtained from the 
eigenvalue problem of a dynamical matrix D(q), 


E D a p(KK,q)W 0 (K,cij ) =ul j W a (K,qj), (6) 

k' /3 


with 

D aP ( KK \ q ) = -^= v$^(o 

y/m K m K > y 

( 7 ) 

Using Eq. ([5]) the hamiltonian Jt? is easily expressed us¬ 
ing creation and annihilation operators. We obtain the 
harmonic part as a sum of harmonic oscillators, 

flUIqj (— + Uqj 'j , (8) 

qj k / 


and the third-order potential as a sum of 3-phonon colli¬ 
sions 

<^3 = y + at A )(fix' + + at*.')- 

AA'A" 

( 9 ) 
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Here the phonon modes (q, j) and (—q, j) have been ab- the strength of interaction between the three phonon A, A' 
breviated by A and —A, respectively. <f> AA / A // represent and A' involved in the scattering. c E > aa'A" is explicitly 

given by 
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where A(q + q'+ q") = 1 when q + q' + q" is a reciprocal 
lattice vector, and is zero otherwise. In Eqs. |[7|) and 
(101 we have used the translational invariance condition 
to simplify the more symmetric equations which appears 
in Ref. [lS- 


The imaginary part of the self energy r A (w) can now 
be computed up to second order in using many body 
perturbation theory. It takes a form analogous to the 
Fermi’s golden rulep^ 


187T 


T A (w) = -p- 55 |$— AA'A"j {(«A' + n \" + 1 )<$(<*> - WA' - WA") + {n\> - n\») [<5(w + CJy - ^A") - S(uj - U\' +WA")]} , 


A'A" 


( 11 ) 


where n A is the phonon occupation number at the equi¬ 
librium, 


n\ = 


exp(hoj\/k B T) - 1' 


( 12 ) 


The double Brillouin zone summation in Eq. (11) is re¬ 
duced to a single summation due to A(—q + q r + q") in 
the calculation. 

2T A (u; A ) corresponds to the phonon linewidth of the 
phonon mode A and its reciprocal is known as the phonon 
lifetime,^ 


1 

2r A ( WA )‘ 


(13) 


B. Lattice thermal conductivity 

When the LBTE is solved under the SMRT method, 
the lattice thermal conductivity tensor can be written in 
a closed formal 

K = 55 CaVa ® v at a mrt , (14) 

where Vo is the volume of a unit cell, v A and r | MRT are 
the group velocity and SMRT of the phonon mode A, 
respectively. C\ is the mode dependent heat capacity 


defined as 


Cx 


k B 


( V exp(hu}\/k B T) 

V k B T J [exp (huj x /k B T) - l] 2 ' 


(15) 


The group velocity can be obtained directly from the 
eigenvalue equation ([6]), 


v a (X) = 


duJx 

dq a 


(16) 




k,k' (3^ 


dq a 


(17) 


In this study it is further assumed that the phonon re¬ 
laxation time is given by the phonon lifetime of Eq. 


SMRT _ ^ 
T X = T X- 


(18) 


This is an approximation. Therefore we have also com¬ 
puted the lattice thermal conductivity from a full solution 
of the LBTE using the method of Ref. CD which is summa¬ 
rized in Appendix [E] Those calculations are much more 
expensive than the closed expression (14). The results 
presented in the Appendix [i] show however that they give 
thermal conductivities which are numerically quite simi¬ 
lar to those with the SMRT approximation for the class 
of compounds we consider. 
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C. Computational details 


Second- and third-order force constants were obtained 
by the supercell approach with finite atomic displace¬ 
ments of 0.03 A. The details of the method are written in 
Appendix [A] and in Ref. HE Non-analytical term correc¬ 
tion 22 was applied to the second-order force constants to 
take into account the long range Coulomb forces present 
in ionic crystals. The technical details to incorporate the 
non-analytical correction are presented in Appendix [B] 
and in Ref. [23] 

For the first-principles calculations, we employed the 
plane-wave basis projector augmented wave method 24 
within the framework of density functional theory as im¬ 
plemented in the VASP codePw^The generalized gradi¬ 
entapproximation (GGA) of Perdew, Burke, and Ernzer- 
hoPSl was U sed for the exchange correlation potential. A 
plane-wave energy cutoff of 500 eV was employed. Recip¬ 
rocal spaces of conventional unit cells of the zincblende- 
type structures and unit cells of the wurtzite-type struc¬ 
tures were sampled by 4 x 4 x 4 and 6x6x4 meshes, 
respectively. The sampling meshes were shifted from 
T-point by (1/2,1/2,1/2) and (0,0,1/2) of neighboring 
mesh grid points along axes, respectively. To obtain 
atomic forces, the total energies were minimized until the 
energy convergences to be less than 10 -8 eV. The lattice 
parameters of the compounds were optimized under the 
calculations at zero pressure condition. The results are 
summarized in Appendix |F| 

The second-order force constants were computed us¬ 
ing 4x4x4 supercells of the conventional unit cell 
(512 atoms) for the zincblende-structures and 5x5x3 
supercells for the wurtzite-type structures (300 atoms). 
Reciprocal spaces of the supercells of the zincblende- 
and wurtzite-type structures were sampled by 1 x 1 x 1 
of T-centered meshes and 1x1x2 meshes shifted by 
(0,0,1/2) with respect to the grid spacing, respectively. 
For zincblende-type GaSb, InN, InAs, and InSb, the sam¬ 
pling meshes were shifted by (1/2,1/2,1/2) to avoid sam¬ 
pling their T-points since their valence and conduction 
bands slightly touch around the T-points due to the un¬ 
derestimation of band gaps by the GGA. Static dielectric 
constant tensors and Born effective charge tensors were 
obtained from density function al per turbation theory as 
implemented in the VASP code.® 9 30 These tensors were 
symmetrized and a sum rule was applied to the Born 
effective charge tensors following Ref. HE The phonon 
band structures and densities of states obtained from 
those harmonic calculation are presented in Appendix |G| 
for all the compounds considered in this study. 

For the third-order force constants, 2x2x2 super¬ 
cells were built for the zincblende-structures (64 atoms) 
and 3x3x2 for the wurtzite-type structures (72 atoms). 
Reciprocal spaces of the supercells of the zincblende- and 
wurtzite-type structures were sampled by 2 x 2 x 2 meshes 
shifted by (1/2,1/2,1/2) and (0,0,1/2) with respect to 
the grid spacing, respectively. These anharmonic calcula¬ 
tions allow to obtain the lifetimes according to Eq. ©■ 


The lifetimes for the phonon states of the band structures 
are shown in Appendix [G] for all the compounds consid¬ 
ered in this study. To obtain those results Brillouin zones 
were sampled using the 96 x 96 x 96 and 80 x 80 x 48 
sampling meshes for the reciprocal spaces of the primi¬ 
tive cells of the zincblende- and wurtzite-type structures, 
respectively. 

Finally, to compute lattice thermal conductivities, the 
reciprocal spaces of the primitive cells of the zincblende- 
and wurtzite-type structures were sampled using the 
19 x 19 x 19 and 19 x 19 x 15 meshes, respectively. In 
Brillouin zone integrations, crystal symmetry was used 
to improve numerical accuracy and to reduce the com¬ 
putational cost. A tetrahedron method was employed 
in the calculation of the imaginary parts of the self 
energies J22HM1 The technical details on the Brillouin zone 
integration are summarized in Appendix 0 The com¬ 
puted lattice thermal conductivities are presented in Ap¬ 
pendix |T] 


III. RESULTS AND DISCUSSION 
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FIG. 3. (color online) Positions of elements in CuA', ZnA, 
and GaA on the periodic table. 


In this section, we attempt to capture the features how 
phonons interact in Brillouin zone schematically by lim¬ 
iting our discussion to selected compounds. In Fig. [T] we 
see a trend that the compounds with larger unit-cell vol¬ 
umes have smaller lattice thermal conductivities. How¬ 
ever the copper compounds, CuCl, CuBr, Cul (CuA), 
clearly ignore this trend. Therefore, to understand this 
behavior, as shown in Fig. [3j we focus on CuA' and their 
neighbors in the periodic table, ZnS, ZnSe, ZnTe (ZnA) 
and GaP, GaAs, GaAs (GaA), and compare them to 
each others to find the characteristics they share. Fig¬ 
ures similar to those found in this section are presented 
in Appendix [G] and [H] for all the zincblende-type com¬ 
pounds considered in this study. 

The phonon band structures between T and L points 
for CuA, ZnA, and GaA are shown in Fig. [4] from the 
left column to the right column, respectively. The pan¬ 
els are aligned in relation to the periodic table in Fig. [3] 
The compounds at the lower panels have heavier anion 
atoms, therefore as it is expected from Eqs. © and (JtJ), 


















5 



FIG. 4. (color online) Phonon linewidth distributions of CuX, ZnX, and Ga.Y of the zincblende-type structure which are 
drawn on the phonon band structures between the F and L points. The scales of the phonon linewidths and phonon frequencies 
are shared in THz, but to make the phonon linewidths visible, their magnitudes are adjusted by scale factors (x0.25, xl, x5, 
and xlO) that are shown close to their compound names. Four phonon branches correspond to TA (blue), LA (green), TO 
(red), and LO (pale blue) modes in order of increasing phonon frequency near F-point, respectively. 


phonon frequencies become systematically lower in the 
lower panels. Looking at the panels from left to right, 
phonon frequencies increase gradually. Therefore, except 
for CuCl and CuBr, the compounds show rather similar 
band structures, at different scales. This means that in 
Eq. (11) the three phonon collisions and decay processes 
allowed by the conservation laws for energy and momen¬ 
tum will be about the same. Then, if the phonon-phonon 
interaction strength can be approximated by a function 
with a simple form, we can expect to analyze those com¬ 
pounds systematically using a small number of parame¬ 
ters. 


CuCl and CuBr reveal an interesting behavior. For all 
compounds there are four phonon branches between T 
and L points. They are labeled as TA (transverse acous¬ 
tic), LA (longitudinal acoustic), TO (transverse optical), 
and LO (longitudinal optical) branches in ascending or¬ 
der along vertical axis near T point. Usually, acoustic 
branches have larger dispersion and group velocity than 
optical branches and the phonon mode frequency of an 
optical branch is weakly dependent on wave vector. How¬ 


ever CuCl shows relatively large dispersions in the optical 
branches. In non-metallic crystals, atomic displacements 
of LO-like phonon modes near T point induce a macro¬ 
scopic electric field. As a result, the phonon mode fre¬ 
quency of an LO branch is lifted up and LO and TO 
branches split near T point. Therefore phonon mode fre¬ 
quency of an LO branch near T is higher than that near 
Brillouin zone boundary in the most compounds. Clear 
LO-TO splitting is observed in CuCl and CuBr, however 
the slopes of the LO branches are inverted. 


Inverse phonon lifetimes, or linewidths, given in 
Eq. (131 are drawn on the phonon band structures as 
the vertical widths of the shaded areas in Fig. [4} Larger 
linewidth indicates larger phonon scattering and shorter 
phonon lifetime. The physical unit of THz is shared 
with phonon mode frequency. To well visualize distri¬ 
butions of the linewidths, the scales of the linewidths 
are adjusted with multiplying factors. It has to be no¬ 
ticed that phonon lifetimes can be different in logarith¬ 
mic scale among compounds, which is presumed from 
the fact that, as shown Fig. [l] lattice thermal conduc- 
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tivities also range in logarithmic scale. It is found that 
linewidths vary non-smoothly with respect to wave vec¬ 
tor. Partly, this come from the conservation laws for 
energy and momentum which select regions where more 
scattering mechanism are possible. As mentioned before, 
this would also explain the similarities in lifetime distri¬ 
bution for compounds with similar band structures. This 
will be described in more details. 


A. (q, uj) map of self energy 


To find common features, linewidth distributions in 
Fig- [4] are observed: linewidths of the TA branches are 
small compared with the other branches, knot-like distri¬ 
butions are found on the LA branches, and, except for 
CuCl, the TO branch linewidths are roughly constant. 
To understand how this variety of the linewidth distri¬ 
butions is obtained, we investigated the frequency and 
wave vector dependencies of the imaginary parts of the 
self energies. The result of these calculations for all the 
zincblende-type compounds considered in this study are 
collected in Appendix [H| 

Taking ZnTe as an example, we draw (q, uj) maps of 
the imaginary parts of the self energies for the TA, LA, 
TO, and LO branches. As shown in Fig. [5j they are 
drawn between T and L points for uj < 6 THz. The 
four branches of phonon band structure are superimposed 
on the corresponding (q, uj) maps as solid curves. The 
linewidth distribution of ZnTe in Fig. [4] is reproduced by 
collecting 2 Fa(w) values on the solid curves where uj = 
oj\, i.e. 2 Fa(wa)- As a general tendency, it is found that 
the imaginary part of the self energy is weakly dependent 
on wave vector, however it varies vividly along frequency. 

To see their frequency dependency in detail, the pro¬ 
files of the (q, uj) maps at L point are plotted in Fig. [b] 
(a). They are drawn up to w ~ 12 THz since 2 Ta(u;) can 
be nonzero up to twice the maximum phonon mode fre¬ 


quency due to the term S(uj—ujy—ujy') in Eq. (11). Note 


that 10 times different scales are used to plot 2 Ta(w) be¬ 
low and above 6 THz in Fig. 6 ] (a). The profiles in the 
lower frequency range are much smaller than those in the 
high frequency range. For each profile, the strongest peak 
appears at uj ~ 10 THz except for the TA branch. The 
peaks at cj ~ 3.5 THz are relatively large in the frequency 
range where the phonon modes exist, uj < 6 THz. The 
phonon band structure of the LA branch crosses this peak 
frequency as can be seen clearly in Fig. [5] This causes 
the knot-like shape of the LA branch linewidth. 

These peaks of the profiles are mainly brought about 
by the three phonon selection rules. In Eq. © , the 
energy conservation among three phonons is represented 
by the S functions and the conservation of momentum 
appears in Eq. (10 1 with the A function. To understand 


the role of the selection rules to imaginary part of self 
energy, the joint density of states (JDOS) is employed. 
This function is temperature independent. There are two 
classes of three phonon scattering events. The collision 
processes, which correspond to that S(uj — u>\ 1 + wa 2 ) 
or -Fc^Ai i n Eq. (Ill, and the decay processes 

which correspond to 6(uj — uj\ 1 — uj\ 2 ). In Ref.(201 they are 
respectively named class 1 and class 2 events, and we will 
keep this terminology in the following to label different 
quantities. For example the JDOS Z? 2 (q, uj) is defined by 
the sum of the JDOS of the two classes and 

-C4 2) (q,w) as 


A>( q,w) = -^(q.w) + D^ ; (q,w) 


( 2 )/ 


(19) 


where 


1 


D 2 (q> w) = — A(-q + q' + q") 


A'A" 

x [<5(w + ujy - W\") + S(uj - uj\> + W A ")] - 
for class 1 events and 

d ( 2 ] (q, w ) = 4 A(-q + q' + q") 


A'A" 


X S(uJ - UJy - UJ A"), 


for class 2 events. Z) 2 (q,u;), D^\q,uj), and Z) 2 “\q,u;) 
are shown in Fig. [ 6 ] (b) for point L. The peak positions 
of the profiles coincide with those of the imaginary parts 
of the self energies. Moreover uj) and D^\q,uj) 

cover the low and high frequency ranges, respectively. 
Therefore it is expected that the class 1 (2) event tends 
to contribute to the phonon lifetimes of the acoustic (op¬ 
tical) branches. 

Comparing Figs, [b] (a) and (b), relative intensities of 
the profiles are clearly different since phonon-phonon in¬ 
teraction strength and phonon occupation numbers are 
ignored in the JDOS. If we assume that the scattering 
of the phonon A with A' and A" can be averaged over 
X' and X" then the imaginary part of self energy can be 
approximated as 


r qjM = -p-P_qjA r 2 (q,w) (20) 

with 

^ = <*b*5>”'A"l 2 . ( 21 ) 

V a ) yy, 

N 2 (q, uj) = 1V 2 (1) (q, uj) + 1V 2 (2) (q, uj) (22) 

and 







7 



FIG. 5. (color online) Imaginary parts of self energies of ZnTe as functions of wave vector between F and L points and frequency 
at 300 K for TA, LA, TO, and LO branches. Phonon dispersion curves of the corresponding branches are superimposed. 


iV 2 (1) (q, W ) 

ivfW) 


^2 A(-q + q' + q")(«A' - n\»)[6(u} + - u)\») - 8{u - u)\> + wa")] 

v A'A" 

^2 A(-q + q' + q")(^A' + «A" + 1 )<5(w - w\> - uj\»), 

V A'A" 


(23) 

(24) 


where n a in Eq. (21) is the number of atoms in the prim¬ 
itive cell. 


The computational cost of Ta(w) may be much smaller 
than that of T>,(u;) if we can avoid the computation of 
third-order force constants using a good approximation. 
In Eq. (|20b, Pqj represents the averaged interaction expe¬ 
rienced by phonon q j while propagating in the medium. 
7V 2 (q, w) is a weighted JDOS (w-JDOS) which measures 
the number of those interactions, N^\ q,w) of them be¬ 
ing of class 1, and ( q, w) being of class 2. 


In Fig. [6] (c), iV 2 (q, w) at L point is presented at 
300 K. Below w ~ 2 THz, the divergence behavior in 
D < 2 \ q,w) has been canceled by the occupation functions 
in N^ (q, ui) but the intensity of the peak at w ~ 10 THz 
is still comparable to those of the other peaks, which is a 
noticeable difference with 2T^(w). However below uj ~ 6 
THz, where the phonon modes do exist, N 2 ( q, uj) reason¬ 
ably reproduces the shapes of the profiles of 2 Ta(w). 

One consequence is that IV 2 (q, ui) may be used for 
a screening of the lattice thermal conductivity of com¬ 


pounds. The (q, ui) map of iV 2 (q, ui) for ZnTe is shown 
in Fig. [TJ The phonon band structure is superimposed on 
this (q, ui) map. The intensity distribution is reasonably 
similar to those of the imaginary parts of the self energies 
in Fig. [5] Therefore low (high) lattice thermal conductiv¬ 
ity is expected if phonon band structure crosses (avoids) 
the high intensity region with larger group velocity. Such 
a discussion for the lattice thermal conductivity will be 
made elsewhere. 


B. Averaged linewidth 

To sort the linewidth distributions quantitatively, we 
define an averaged linewidth over the Brillouin zone, 
2 T ave , as a temporal measure, with 

Fave= 3nJV^ rA(WA) - (25) 

Averaged linewidths of CuA', Z 11 X, and GaAT are listed 
in Table |T] The averaged linewidths increase from right 
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FIG. 6. (color online) Imaginary parts of self energies (2r), 
JDOS ( D 2 ), and w-JDOS (A 2 ) are drawn for ZnTe as a func¬ 
tion of phonon frequency at L point, (a) Imaginary parts 
of self energies 2r of TA, LA, TO, and LO branches at 300 
K are depicted by dashed, dashed-dotted, solid, and dotted 
curves, respectively. Note that the scale for the vertical axis is 
changed below and above 6 THz. (b) JDOS D ^ and D 2 2 '* are 
drawn by the solid and dotted curves, respectively, and the 
sum D 2 = is depicted as the shaded region, (c) w- 

JDOS V 2 (1) and JV 2 2 ^ at 300 K are presented by the solid and 
dotted curves, respectively, and the sum N 2 = + A r 2 1 ' is 

depicted as the shaded region. 


to left and from lower to upper in Table [T] Along the 
compounds with fixed cations, relatively large differences 
of the averaged linewidths are found between CuBr-Cul 
and ZnS-ZnSe. We may expect a qualitative change 
of phonon collisions at these points. The averaged 
linewidths of the 33 zincblende-type compounds are sum¬ 
marized in Fig. [8] (filled circles). In between BP-BAs, 
AlAs-AlSb, InP-InAs, and CdS-CdSe, large differences of 
the averaged linewidths are also found. These qualitative 
changes are also observed in the linewidth distributions, 
which are shown in Appendix [Gj 


TABLE I. Averaged linewidths 2T ave of CuA, ZnA, and 
GaA in THz. See Eq. (25 I for the definition of r ave . 


CuCl 

1.3 x 10° 

ZnS 

1.2 x 1CT 1 

GaP 

4.1 x 10“ 2 

CuBr 

4.9 x 10 _1 

ZnSe 

3.3 x 1CT 2 

GaAs 

2.8 x 10“ 2 

Cul 

3.8 x 10“ 2 

ZnTe 

2.5 x 1CT 2 

GaSb 

1.6 x 10“ 2 



FIG. 7. (color online) w-JDOS (IV 2 ) of ZnTe as functions of 
wave vector between T-L and phonon frequency at 300 K. The 
solid curves show phonon mode frequencies of the TA, LA, 
TO, and LO branches in ascending order along the vertical 
axis. 


The averaged linewidths obtained with an averaged 
phonon-phonon interaction, 2T ave are also shown in 
Fig. [i] (open circles). 2r ave tends to overestimate 2r ave , 
however the degree of the overestimation is not system¬ 
atic among the 33 zincblende-type compounds. This sug¬ 
gests that phonon-phonon interaction strength is as com¬ 
plicated as three phonon selection rule and their detailed 
combination is decisively important to determine phonon 
linewidth quantitatively. 


IV. SUMMARY 


The lattice thermal conductivities of the zincblende- 
and wurtzite-type compounds with 33 combinations of 
elements are calculated with RTA and LBTE from first- 
principles anharmonic lattice dynamics calculations. The 
zincblende- and wurtzite-type compounds with the same 
compositions show similar lattice thermal conductivities, 
phonon densities of states, and phonon linewidth distri¬ 
butions though their stacking sequences are different. Fo¬ 
cusing on the zincblende-type compounds and selecting 
CuA, ZnA, and GaA, the phonon linewidth distribu¬ 
tions are investigated in detail. It is found that their 
phonon linewidths are distributed non-smoothly in the 
Brillouin zones. This is explained by the weak wave 
vector dependence and the strong frequency dependence 
of the imaginary parts of the self energies. It is shown 
that detailed combination of phonon-phonon interaction 
strength and three phonon selection rule is critically 
important to determine imaginary part of self energy. 
Therefore it is considered difficult to predict phonon life¬ 
time quantitatively without a full calculation using an¬ 
harmonic force constants. However we show a possibil¬ 
ity that w-JDOS may be used for a rough estimate of 
the lattice thermal conductivity and thereby materials 
screening. 
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FIG. 8. (color online) Averaged linewidths 2r ave (filled circle) and 2r ave (open circle) of the 33 zincblende-type compounds. 
The solid and dashed lines are eye guide to group compounds. 
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Appendix A: Supercell approach and finite 
displacement approach 


We employed the supercell approach and the finite dis¬ 
placement approach to calculate force constants. The 
harmonic and cubic force constants are given by 


$ q/ 3 (Zk,ZV) 


d 2 $ 

du a (ln)dup(l' k') ’ 


(Al) 


and 


$ a/ 3 7 (lK, I'k', l" k") 


<9 3 $ 

du a {ln)dup(V K')du-y(l" k") 1 

(A2) 


respectively. A set of finite atomic displacements is in¬ 
troduced to supercells to calculate these force constants. 
Forces acting on atoms in a supercell are obtained from 
the first-principles calculations and the force constants 
are calculated from the finite difference method as fol¬ 
lows. The harmonic force constants may be approxi¬ 
mately obtained by 


$ a p(I k, I'k') 


Fp[l'n'; u(Zk)] 

U a (lK) 


(A3) 


where F[1'k'; u(Zre)] is defined as the atomic force mea¬ 
sured at r (I'k 1 ) under an atomic displacement u (Ik) in 
a supercell. The cubic force constants may be similarly 
obtained by 


c E > q k', l" k") 


F 1 [1”k"-,u{Ik),u(1'k')\ 
u a (ln)up(l' k') 


where F[l" k"\ u (Ik), u(Z / «; / )] is defined as the atomic force 
measured at r (l" k") under a pair of atomic displacements 
u (Ik) and u (I'k 1 ) in a supercell. The sets of Eqs. (A31 and 
(A4) forms systems of linear equations, respectively, and 
they were solved using the Moore-Penrose pseudoinverse. 
The computational details are found in Ref. [21, 

Configurations of atomic displacements in supercells 
have to be chosen to solve the linear equations. The 
number of the atomic displacement configurations can be 
huge, however it may be reduced utilizing crystal sym¬ 
metry. Atoms that are independent with respect to space 
group symmetry are selected to be displaced in supercells. 
The number of displacement directions is reduced using 
site symmetry at the displaced atom. For the case of 
cubic force constants, the second atom that is displaced 
is determined from the symmetry lowered by the first 
atomic displacement. To gain numerical stability, the 
opposite direction of a chosen atomic displacement was 
introduced if they are symmetrically inequivalent. The 
crystal symmetry is again used to recover full elements 
of force constants from the reduced set of supercell force 
calculations. Let a site-symmetry operation S at the dis¬ 
placed atom r (Ik) send an atomic point to the other point 
as r {I'k') -A r (I'k'). With this site-symmetry operation 


and from Eq. (A3), the first atomic displacement, forces 
on atoms in the supercell, and harmonic force constants 
are related by 


S' q , 7 u 7 (Zk)$ q , / 3(Ik, I'k') ~ Sp-yF^[l'k'; u(Zk)]. 

ol '7 7 


(A5) 

Similarly for the cubic force constants, Eq. (A4) is rewrit¬ 
ten as 


l' k' , l" k") ~ 
c*e 

Sp C S^A^ Cr] [l'k', l"k": u (i/s)], (A6) 
Cv 

where 


A4> / 3 7 [Z'k', 1"k"; u(Zre)] = 

F 7 [1"k"', u (Ik), u(Z / k')] 
up (l' k') 


^ 0j (1'k',1"k"). (A7) 


Appendix B: Non-analytical term correction 


Macroscopic electric field is generated by atomic dis¬ 
placements of phonon modes. At long-wavelength limit, 
this effect is incorporat ed as th e non-analytical part of 
harmonic force constant d 22 * 35 * 36 ! given by 




47 re 2 E 7 < h Z K,~ta Ey < h' Z K' n >p 
Vo E 77 '?7 e “'97' ." 


where e is the elemental charge, Z* is the Born effec¬ 
tive charge tensor for the atom k, and e°° is the high 
frequency static dielectric constant tensor. The non- 
analytical term is q dependent but only on its direction in 
reciprocal space and its application is limited at q A- 0. 
In the supercell approach, the macroscopic electric field 
effect at general q-points may be included into the dy¬ 
namical matrix using a Fourier interpolation scheme pre¬ 
sented in Ref. 12.41 The summation in Eq. 0 is made 
over lattice points in a supercell and the non-analytical 
term is added to the harmonic supercell force constants 
divided by the number of unit cells in the supercell, N c , 
as 


$^ c (^,ZV,q) = $ a p(lK,r K ')+$™ f) (KK',t l )/N c , (B2) 


which is used instead of normal & a p(lK,l' k') to create 
dynamical matrices in Eq. ([7| . 

In Eq. (16), first derivative of dynamical matrix re¬ 
quires first derivative of the non-analytical term that is 
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easily obtained as 
a^O^.q) _ 47re 2 

dqc, Vo 

^kX<x E 7 ' 97'-^k',7'/3 + ^V,C 0 E 7 ( h^K,ja 

E 77 ' 

E 7 Q-y^K,,ja E 7 ' Qj’^k’^'P E 7 " 9 7 "( e ^" + e y'f) 

(E 77 ' <h£™'<h') 

(B3) 

Hence the first derivative of the dynamical matrix with 
the non-analytical term correction becomes 

dD^ c (KK',q) 1 1 

dqc V m K m ' K N ° 

E f 1 C«',q) . r ,x o 

I Q P v zq-[r(Z k ) —r(Z/c)J , 

,r ^ 

i$^ c (^,ZV,q) [r c (ZV) -r c ^)]e iq '[ r(zV) - r(iK) ]|. 

(B4) 


Appendix C: Brillouin zone integration 


Computations of phonon lifetimes were discretized us¬ 
ing sampling meshes in Brillouin zones. The calculation 
order of a phonon lifetime in Eq. (JTTj) is 0(n 6 ) with re¬ 
spect to the number of mesh points along an axis. This is 
reduced using translation symmetry of crystal to 0(n 3 ) 
by A(q + q' + q"). In Eq. (11), it is supposed that the 
leftmost index of the phonon-phonon interaction strength 
(Eq. 10) is fixed and the rest of the indices can be inter¬ 
changed as 4>_aa'A" = < I > _aa"A' • Therefore the calcula¬ 
tion order becomes 0(n 3 / 2). A set of q-point triplets, 
{(q, q', q")| q, q', q" € Brillouin zone} was thus reduced. 
The number of the q-point triplets was reduced further 
considering crystallographic point group symmetry and 
time-reversal symmetry. 

W e p erformed the Brillouin zone in tegrat ion in 
Eq. ml using a linear tetrahedron method ES K a S Smear¬ 
ing method is another popular method that is easier to 
implement in computer programs. In a smearing method, 
the delta functions are simply replaced by the Gaus¬ 
sian or Lorentzian functions that are normalized so that 
those integrals become 1. The advantage of the tetrahe¬ 
dron method over the smearing method is that it can be 
parameter-free. The straightforward application of the 
tetrahedron method to Eq. (11) requires the translational 


invariance condition where Eq. (|ll ) is reduced to a single 
summation over Brillouin zone)®By sorting the order of 
integration of contributions from vertices of tetrahedra, 
each delta function is simply replaced by respective inte¬ 
gration weight that is pre-calculated from phonon mode 


frequencies on a sampling mesh. 33 Therefore after com¬ 
puting the integration weights, the tetrahedron method 
and smearing method can be treated in a computer pro¬ 
gram on equal footing. Generally the smearing method 
with an appropriate smearing parameter underestimates 
phonon lifetime, and thus lattice thermal conductivity as 
well. In both methods, the integrations are necessary to 
be performed over the reduced set of q-points or triplets. 


The Brillouin zone integration of Eq. (14) reduces to 


the integration over an irreducible part of the Brillouin 
zone (IBZ): 

Ka>3 = NVo E 

u qelBZ j 

-EE R a ^RpCV 7 (qj)vc{qj), (Cl) 

R 7C 


where R runs the point group operations and p q is the 
order of point-group of the q-vector. 


Appendix D: Isotope scattering 


The scattering rate for a phonon mode by randomly 
distributed isotopes is given by the second-order pertur¬ 
bation theory ns 37 


1 


w 


A' 


WA') 




2 , (Dl) 


where g K is the mass variance parameter, g K is defined 

by 

9k = E] M 1 ~ m iK./in K ) 2 , (D2) 


where fi and rrii^ are the mole fraction and relative 
atomic mass of the *th isotope, respectively. fh K is the 
average mass given by fh K = fi m iK- When both 
of phonon-phonon and phonon-isotope scattering effects 
were considered in the calculation, we employed the scat¬ 
tering rate given as the sum of them: 


11 1 

r} otal t x + t} so (wa) 


(D3) 


Appendix E: Direct solution of linearized phonon 
Boltzmann equation 

Recently the method of the direct solution of the lin¬ 
earized phonon Boltzmann equation was developed by 
ChaputFb In this method, lattice thermal conductivity is 
given by 

_ ^ _ U\Vg(X) U\’Vp(\') , n ^i. 

P 4kBT2NV ° AA' sinh sinh (jg^) 

(El) 
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where D""' 1 denotes the Moore-Penrose inverse of the col¬ 
lision matrix Q that is given by 


^AA' 



7r 

n? 


EI^AA'A" 

A" 



X [5(u>x-uj' x -u)”) 

+ 5 (wa + u x — Uj'x) + S(uj\ — U! x + U>x) ]. 


(E2) 


To include the phonon -isot ope scattering effect in Ap¬ 
pen dix h| t\ in Eq. (E2| was replaced by rj; otal of 
Eq. ( fmF 


Appendix F: Lattice parameters 

Since phonon lifetime is sensitive to lattice parame¬ 
ters, the lattice parameters and volumes of the (conven¬ 
tional) unit cells of the zincblende- and wurtzite-type 
compounds with 33 combinations of chemical elements 
used in the present study are listed in Table |TT] 


Appendix G: Phonon band structures 

Figures from [9] to |20| show the phonon band structures 
and distributions of phonon linewidths (left panels) and 
phonon densities of states (DOS) (right panels) for the 
zincblende- and wurtzite-type compounds with 33 com¬ 
binations of elements used in the present study. The 
reciprocal path of T-L of the zincblende-type structure is 
around twice longer than that of T-A of the wurtzite-type 
structure. The dotted curves in the left panel show the 
phonon band structure of the zincblende-type compound 
folded at the middle point between L and T points to 
emphasize the similarity between the phonon band struc¬ 
tures of the zincblende- and wurtzite-type compounds. 
The phonon DOS shaded and drawn by solid curves are 
those of the zincblende- and wurtzite-type compounds, 
respectively. 


TABLE II. Lattice parameters (A) and volumes (A 3 ) of con¬ 
ventional unit cells of the zincblende- and wurtzite-type struc¬ 
tures used in the present study. 


Zincblende-type Wurtzite-type 



a 

Volume 

a 

c 

volume 

BN 

3.63 

47.6 

2.56 

4.22 

23.9 

BP 

4.55 

94.0 

3.20 

5.30 

47.0 

BAs 

4.81 

111.1 

3.38 

5.60 

55.5 

AIN 

4.40 

85.2 

3.13 

5.02 

42.5 

A1P 

5.51 

166.9 

3.89 

6.38 

83.5 

AlAs 

5.73 

187.8 

4.04 

6.64 

93.9 

AlSb 

6.23 

242.1 

4.40 

7.23 

121.0 

GaN 

4.55 

94.1 

3.22 

5.24 

47.0 

GaP 

5.51 

166.9 

3.88 

6.39 

83.5 

GaAs 

5.75 

189.9 

4.05 

6.67 

95.0 

GaSb 

6.22 

240.2 

4.38 

7.21 

120.1 

InN 

5.05 

128.9 

3.58 

5.79 

64.4 

InP 

5.96 

211.5 

4.21 

6.90 

105.8 

InAs 

6.19 

236.8 

4.37 

7.17 

118.5 

InSb 

6.63 

291.9 

4.68 

7.69 

146.0 

BeO 

3.83 

56.1 

2.71 

4.41 

28.1 

BeS 

4.88 

115.9 

3.44 

5.67 

58.1 

BeSe 

5.18 

139.2 

3.65 

6.02 

69.7 

BeTe 

5.67 

182.1 

4.00 

6.59 

91.1 

MgTe 

6.51 

276.3 

4.61 

7.50 

138.2 

ZnO 

4.63 

99.2 

3.29 

5.30 

49.6 

ZnS 

5.45 

161.7 

3.85 

6.31 

80.9 

ZnSe 

5.74 

189.2 

4.05 

6.65 

94.7 

ZnTe 

6.18 

236.5 

4.36 

7.18 

118.4 

CdS 

5.94 

209.3 

4.20 

6.84 

104.7 

CdSe 

6.21 

239.2 

4.39 

7.17 

119.7 

CdTe 

6.63 

290.8 

4.68 

7.67 

145.6 

CuCl 

5.43 

159.7 

3.82 

6.32 

79.9 

CuBr 

5.73 

187.7 

4.04 

6.65 

93.9 

Cul 

6.08 

224.8 

4.29 

7.06 

112.4 

CuH 

4.02 

64.9 

2.87 

4.56 

32.5 

Agl 

6.61 

288.8 

4.67 

7.64 

144.3 

SiC 

4.38 

84.0 

3.09 

5.07 

42.0 












13 



FIG. 9. (color online) Linewidth distributions and phonon 
densities of states of the zincblende- and wurtzite-type BN, 
BP, and BAs. The linewidths are drawn by 50 times magni¬ 
fied. 



FIG. 10. (color online) Linewidth distributions and phonon 
densities of states of the zincblende- and wurtzite-type AIN, 
A1P, AlAs, and AlSb. The linewidths are drawn by 20 times 
magnified. 
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FIG. 11. (color online) Linewidth distributions and phonon 
densities of states of the zincblende- and wurtzite-type GaN, 
GaP, GaAs, and GaSb. The linewidths are drawn by 10 times 
magnified. 


FIG. 12. (color online) Linewidth distributions and phonon 
densities of states of the zincblende- and wurtzite-type InN, 
InP, InAs, and InSb. The linewidths are drawn by 10 times 
magnified. 
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Zincblende Wurtzite DOS 



FIG. 13. (color online) Linewidth distributions and phonon 
densities of states of the zincblende- and wurtzite-type BeO, 
BeS, BeSe, and BeTe. The linewidths are drawn by 20 times 
magnified. 


FIG. 15. (color online) Linewidth distributions and phonon 
densities of states of the zincblende- and wurtzite-type ZnO, 
ZnS, ZnSe, and ZnTe. The linewidths are drawn by 5 times 
magnified. 



FIG. 14. (color online) Linewidth distributions and phonon 
densities of states of the zincblende- and wurtzite-type MgTe. 
The linewidths are drawn by 5 times magnified. 
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FIG. 16. (color online) Linewidth distributions and phonon 
densities of states of the zinc.blende- and wurtzite-type CdS, 
CdSe, and CdTe. The linewidths are drawn by 5 times mag¬ 
nified. 



FIG. 18. (color online) Linewidth distributions and phonon 
densities of states of the zincblende- and wurtzite-type CuBr, 
Cul, and CuH. 




FIG. 17. (color online) Linewidth distributions and phonon FIG. 19. (color online) Linewidth distributions and phonon 
densities of states of the zincblende- and wurtzite-type GuCl. densities of states of the zincblende- and wurtzite-type Agl. 
The linewidths are drawn by 4 times reduced. 
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FIG. 20. (color online) Linewidth distributions and phonon 
densities of states of the zincblende- and wurtzite-type SiC. 
The linewidths are drawn by 20 times magnified. 


Appendix H: Imaginary parts of self energies 


Figures from 21 to 53 present (q, to) maps of imaginary 


parts of self energies and w-JDOS of the 33 zincblende- 
type compounds used in the present study. Their phonon 
band structures are superimposed on the (q, lj) maps. 
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FIG. 21. (color online) (q, cj) maps of imaginary parts of self 
energies 2 Fa(^) and w-JDOS iV 2 (q, cj) of BN. 
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FIG. 22. (color online) (q, cj) maps of imaginary parts of self 
energies 2Fa(cj) and w-JDOS A^q, cj) of BP. 


FIG. 23. (color online) (q, cj) maps of imaginary parts of self 
energies 2 Fa(cj) and w-JDOS A^q, cj) of BAs. 




























































19 



0.080 

0.064 


0.048 U 
X 

H 

0.032 U 


0.016 

0.000 


0.5 

0.4 

°- 3 S 

H. 

- 0.1 

L-Io.o 

0.54 

0.48 

0.42 

0.36 

0.30 g 

0.24 ~ 

0.18 M 

0.12 

0.06 

0.00 

0.40 

0.32 

0.24 ' S ' 

S 

0.16 u 

- 0.08 
0.00 



» 0.09 
0.08 
0.07 
0.06 
_ 0.05 
- 0.04 
- 0.03 
- 0.02 
- 0.01 
Jo.oo 


£ 





■ °- 20 
I 0.16 

B 0.12 

■ X 

H 

- 0.08 u 

- 0.04 

Jo.oo 



» 0.225 
0.200 
0.175 
0.150 
0.125 x 
- 0.100 ~ 
- 0.075 04 
- 0.050 
- 0.025 

-I 0.000 




FIG. 24. (color online) (q, uj) maps of imaginary parts of self 
energies 2 Fa(^) and w-JDOS iV 2 (q, cj) of AIN. 


FIG. 25. (color online) (q,cj) maps of imaginary parts of self 
energies 2 Fa(^) and w-JDOS iV 2 (q, cj) of A1P. 
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FIG. 26. (color online) (q, cu) maps of imaginary parts of self 
energies 2 Fa(ci;) and w-JDOS A^(q, cj) of AlAs. 


FIG. 27. (color online) (q, uj) maps of imaginary parts of self 
energies 2 Fa(o;) and w-JDOS A^q, co) of AlSb. 
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FIG. 28. (color online) (q, uj) maps of imaginary parts of self 
energies 2 Fa(^) and w-JDOS iV 2 (q, uj) of GaN. 


FIG. 29. (color online) (q, uj) maps of imaginary parts of self 
energies 2 Fa(^) and w-JDOS iV 2 (q, uj) of GaP. 
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FIG. 30. (color online) (q, uj) maps of imaginary parts of self 
energies 2 Fa(cj) and w-JDOS iV 2 (q, oj) of GaAs. 


FIG. 31. (color online) (q, cj) maps of imaginary parts of self 
energies 2 Fa(cj) and w-JDOS A^(q, cj) of GaSb. 
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FIG. 32. (color online) (q, cu) maps of imaginary parts of self 
energies 2 Fa(ci;) and w-JDOS A^(q, co) of InN. 


FIG. 33. (color online) (q, uj) maps of imaginary parts of self 
energies 2 Fa(cj) and w-JDOS A^q, oj) of InP. 
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FIG. 34. (color online) (q, uj) maps of imaginary parts of self 
energies 2 Fa(^) and w-JDOS A^ 2 (q, oj) of InAs. 


FIG. 35. (color online) (q,cj) maps of imaginary parts of self 
energies 2 Fa(^) and w-JDOS A^q, oj) of InSb. 
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FIG. 36. (color online) (q, uj) maps of imaginary parts of self 
energies 2 Fa(^) and w-JDOS A^q, uj) of BeO. 


FIG. 37. (color online) (q, uj) maps of imaginary parts of self 
energies 2 Fa(^) and w-JDOS A^q, uj) of BeS. 


W 2 (THz - ') 2r (THz) 2T (THz) 2T (THz) 





























































26 



0.060 
0.048 
0.036 1 ? 


0.012 

0.000 


(b) 


EH 

H 


(C) 


H 


(d) 


EE 

H 



0.225 
0.200 
0.175 
0.150 
0.125 x 
0.100 ~ 
0.075 ^ 
0.050 
0.025 
0.000 



- 0.3 
0.0 



- 0.3 
0.0 




» 0.054 

0.048 

0.042 

0.036 

0.030 £ 
H 

- 0.024 ~ 
- 0.018 M 
- 0.012 
- 0.006 
-10.000 

» 0.135 
0.120 
0.105 
0.090 
0.075 x 
- 0.060 ~ 
- 0.045 01 
- 0.030 
- 0.015 

J 0.000 



1 0.54 
0.48 
0.42 
0.36 
0.30 £ 
0.24 ^ 
0.18 04 
0.12 
0.06 
0.00 



» 0.45 
0.40 
0.35 
0.30 
0.25 
- 0.20 
- 0.15 
- 0.10 
- 0.05 

Jo.oo 




FIG. 38. (color online) (q, uj) maps of imaginary parts of self 
energies 2 Fa(^) and w-JDOS A^q, uj) of BeSe. 


FIG. 39. (color online) (q, uj) maps of imaginary parts of self 
energies 2 Fa(^) and w-JDOS iV 2 (q, uj) of BeTe. 
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FIG. 40. (color online) (q, uj) maps of imaginary parts of self 
energies 2 Fa(cj) and w-JDOS iV 2 (q, cj) of ZnO. 


FIG. 41. (color online) (q, cj) maps of imaginary parts of self 
energies 2 Fa(cj) and w-JDOS A^ 2 (q, cj) of ZnS. 
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FIG. 42. (color online) (q, uj) maps of imaginary parts of self 
energies 2 Fa(^) and w-JDOS A^q, uj) of ZnSe. 


FIG. 43. (color online) (q, uj) maps of imaginary parts of self 
energies 2 Fa(^) and w-JDOS A^ 2 (q, uj) of ZnTe. 
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FIG. 44. (color online) (q, u) maps of imaginary parts of self 
energies 2r.\(a;) and w-JDOS A r 2 (q. cd) of CdS. 


FIG. 45. (color online) (q, u) maps of imaginary parts of self 
energies 2FA(a;) and w-JDOS A r 2 (q. cj) of CdSe. 
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FIG. 46. (color online) (q, uj) maps of imaginary parts of self 
energies 2 Fa(^) and w-JDOS iV 2 (q, uj) of CdTe. 


FIG. 47. (color online) (q, uj) maps of imaginary parts of self 
energies 2 Fa(^) and w-JDOS iV 2 (q, uj) of MgTe. 
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FIG. 48. (color online) (q, uj) maps of imaginary parts of self 
energies 2 Fa(^) and w-JDOS A^ 2 (q, uj) of CuCl. 


FIG. 49. (color online) (q, uj) maps of imaginary parts of self 
energies 2 Fa(^) and w-JDOS A^ 2 (q, uj) of CuBr. 
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FIG. 50. (color online) (q, oj) maps of imaginary parts of self 
energies 2r_\(m) and w-JDOS ^(q, oj) of Cul. 


FIG. 51. (color online) (q, oj) maps of imaginary parts of self 
energies 2 Fa ( u ;) and w-JDOS ^(q, to) of CuH. 
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FIG. 52. (color online) (q, uj) maps of imaginary parts of self 
energies 2 Fa(^) and w-JDOS A^q, uj) of Agl. 


FIG. 53. (color online) (q, uj) maps of imaginary parts of self 
energies 2 Fa(^) and w-JDOS A^(q, cj) of SiC. 
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Appendix I: Lattice thermal conductivities 


Lattice thermal conductivities of the 33 zincblende- and wur tzite -typ e compounds calculated with RTA and the 
method of the direct solution of LBTE^ are listed in Tables [ill] and |vj Scattering by the isotopes with the natural 
abundance 38 is considered in Table IV Available experimental values in reports are shown in Table IV and comparisons 
between the experiments and calculations are graphically shown in Fig. [2] The methods to include isotope effect and 
the direct solution of LBTE are described in Appendices [D] and [Ej respectively. 


TABLE III. Calculated lattice thermal conductivities of the 33 zincblende- and wurtzite-type compounds at 300 K in W/m-K. 


Zincblende-type 



Wurtzite-type 



RTA 

LBTE 


RTA 


LBTE 


XX 

XX 

XX 


ZZ XX 


zz 


BN 

1.17x 

10 3 

1.80xl0 3 

9.30 

X 

10 2 

9.24 

X 

10 2 

1.23 

X 

10 3 

1.04 

X 

10 3 

BP 

4.27x 

10 2 

4.64xl0 2 

4.18 

X 

10 2 

3.19 

X 

10 2 

4.31 

X 

10 2 

3.28 

X 

10 2 

BAs 

1.59x 

10 3 

5.61xl0 3 

1.44 

X 

10 3 

1.09 

X 

10 3 

2.38 

X 

10 3 

1.21 

X 

10 3 

AIN 

2.12x 

10 2 

2.35x 10 2 

2.24 

X 

10 2 

2.07 

X 

10 2 

2.40 

X 

10 2 

2.12 

X 

10 2 

A1P 

8.30x 

10 1 

8.40X10 1 

7.48 

X 

10 s 

7.14 

X 

10 s 

7.46 

X 

10 s 

7.17 

X 

10 s 

AlAs 

8.56x 

10 1 

8.68X10 1 

7.23 

X 

10 s 

6.96 

X 

10 s 

7.39 

X 

10 s 

7.10 

X 

10 s 

AlSb 

8.92 x 

10 1 

9.48X10 1 

6.18 

X 

10 s 

5.84 

X 

10 s 

6.35 

X 

10 s 

5.82 

X 

10 s 

GaN 

2.31x 

10 2 

3.37X10 2 

2.56 

X 

10 2 

2.38 

X 

10 2 

2.73 

X 

10 2 

2.66 

X 

10 2 

GaP 

9.62x 

10 1 

1.22x 10 2 

1.04 

X 

10 2 

9.08 

X 

10 s 

1.11 

X 

10 2 

9.63 

X 

10 s 

GaAs 

3.18x 

10 1 

3.31X10 1 

2.85 

X 

10 s 

2.61 

X 

10 s 

2.86 

X 

10 s 

2.66 

X 

10 s 

GaSb 

3.30x 

10 1 

3.39X10 1 

2.42 

X 

10 s 

1.97 

X 

10 s 

2.48 

X 

10 s 

2.02 

X 

10 s 

InN 

9.54x 

10 1 

1.08xl0 2 

1.10 

X 

10 2 

1.13 

X 

10 2 

1.18 

X 

10 2 

1.19 

X 

10 2 

InP 

7.51 x 

10 1 

8.71X10 1 

6.87 

X 

10 s 

6.57 

X 

10 s 

7.02 

X 

10 s 

6.93 

X 

10 s 

InAs 

2.05x 

10 1 

2.53X10 1 

1.83 

X 

10 s 

1.71 

X 

10 s 

1.86 

X 

10 s 

1.80 

X 

10 s 

InSb 

1.27x 

10 1 

1.40X10 1 

9.71 

X 

10° 

9.63 

X 

10° 

9.83 

X 

10° 

9.93 

X 

10° 

BeO 

3.58x 

10 2 

4.64x10 s 

2.86 

X 

10 2 

2.98 

X 

10 2 

3.36 

X 

10 2 

3.20 

X 

10 2 

BeS 

1.57x 

10 2 

1.74x10 s 

1.29 

X 

10 2 

1.21 

X 

10 2 

1.37 

X 

10 2 

1.27 

X 

10 2 

BeSe 

3.89x 

10 2 

6.75X10 2 

3.21 

X 

10 2 

3.12 

X 

10 2 

4.08 

X 

10 2 

3.36 

X 

10 2 

BeTe 

2.86x 

10 2 

5.03x10 s 

1.83 

X 

10 2 

1.82 

X 

10 2 

2.36 

X 

10 2 

1.93 

X 

10 2 

MgTe 

1.57x 

10 1 

1.57x10 s 

1.22 

X 

10 s 

1.36 

X 

10 s 

1.24 

X 

10 s 

1.35 

X 

10 s 

ZnO 

5.82x 

10 1 

8.93x10 s 

4.63 

X 

10 s 

6.24 

X 

10 s 

4.66 

X 

10 s 

6.48 

X 

10 s 

ZnS 

3.72 x 

10 1 

4.72x10 s 

3.73 

X 

10 s 

4.03 

X 

10 s 

3.89 

X 

10 s 

4.69 

X 

10 s 

ZnSe 

1.63x 

10 1 

1.72x10 s 

1.50 

X 

10 s 

1.53 

X 

10 s 

1.51 

X 

10 s 

1.57 

X 

10 s 

ZnTe 

1.96x 

10 1 

2.04x10 s 

1.45 

X 

10 s 

1.55 

X 

10 s 

1.45 

X 

10 s 

1.58 

X 

10 s 

CdS 

2.13x 

10 1 

2.30x10 s 

1.79 

X 

10 s 

2.01 

X 

10 s 

1.88 

X 

10 s 

2.05 

X 

10 s 

CdSe 

l.lOx 

10 1 

1.24x10 s 

8.86 

X 

10° 

9.58 

X 

10° 

8.81 

X 

10° 

1.04 

X 

10 s 

CdTe 

6.51x 

10° 

6.95x10° 

5.03 

X 

10° 

5.67 

X 

10° 

5.01 

X 

10° 

5.87 

X 

10° 

CuCl 

1.48 x 

10° 

1.11x10° 

1.01 

X 

10° 

1.75 

X 

10° 

1.04 

X 

10° 

2.02 

X 

10° 

CuBr 

2.55 x 

10° 

2.79x10° 

1.77 

X 

10° 

2.88 

X 

10° 

1.75 

X 

10° 

3.04 

X 

10° 

Cul 

6.53x 

10° 

7.22x10° 

5.89 

X 

10° 

6.94 

X 

10° 

6.07 

X 

10° 

7.32 

X 

10° 

CuH 

1.44 x 

10 1 

1.50x10 s 

8.17 

X 

10° 

6.47 

X 

10° 

8.62 

X 

10° 

8.12 

X 

10° 

Agl 

2.02x 

10° 

2.47x10° 

1.39 

X 

10° 

1.78 

X 

10° 

1.38 

X 

10° 

1.90 

X 

10° 

SiC 

4.15x 

10 2 

4.42xl0 2 

4.05 

X 

10 2 

3.44 

X 

10 2 

4.20 

X 

10 2 

3.48 

X 

10 2 














35 


TABLE IV. Calculated lattice thermal conductivities of the 33 zincblende- and wurtzite-type compounds with isotope effect 
at 300 K in W/m-K. Experimental values at room temperature are also presented. 


Zincblende-type Wurtzite-type 



RTA 

LBTE 

Exp. 

RTA 





LBTE 



Exp. 


XX 

XX 


XX 

zz 


XX 


ZZ 



BN 

6.48xl0 2 

7.26x 10 2 

7.6x 10 2 

* 5.61 x 10 2 

5.42 

X 

10 2 

6.02 

X 

10 2 

5.73 

X 

10 2 

- 

BP 

3.96X10 2 

4.20x 10 2 

3.5x 10 2 

" 3.89 x 10 2 

3.01 

X 

10 2 

3.95 

X 

10 2 

3.07 

X 

10 2 

- 

BAs 

1.23x10 s 

2.94x10 s 

- 

1.15 x 10 s 

9.07 

X 

10 2 

1.73 

X 

10 s 

9.94 

X 

10 2 

- 


AIN 

2.11xl0 2 

2.35x 10 2 

- 

2.24 x 10 2 

2.07 

X 

10 2 

2.39 

X 

10 2 

2.12 

X 

10 2 

3.5xl0 2 [ 

] 

A1P 

8.30X10 1 

8.40X10 1 

9x 10 s F 

7.48 x 10 s 

7.14 

X 

10 s 

7.46 

X 

10 s 

7.17 

X 

10 s 

- 


AlAs 

8.56X10 1 

8.68X10 1 

9.8x10 s 

* 7.23 x 10 s 

6.96 

X 

10 s 

7.39 

X 

10 s 

7.10 

X 

10 s 

- 

AlSb 

7.51X10 1 

7.74x10 s 

5.6x10 s 

5.71 x 10 s 

5.44 

X 

10 s 

5.82 

X 

10 s 

5.40 

X 

10 s 

- 


GaN 

1.54xl0 2 

1.81xl0 2 

- 

1.65 x 10 2 

1.66 

X 

10 2 

1.71 

X 

10 2 

1.72 

X 

10 2 

2.1x10 2 [ 

] 

GaP 

8.55X10 1 

1.04xl0 2 

l.OOxlO 2 

*] 9.19 x 10 s 

8.16 

X 

10 s 

9.65 

X 

10 s 

8.54 

X 

10 s 

- 


GaAs 

3-lOxlO 1 

3.21X10 1 

4.5x 10 s r 

T 2.78 x 10 s 

2.55 

X 

10 s 

2.78 

X 

10 s 

2.59 

X 

10 s 

- 

GaSb 

3.18X10 1 

3.25x10 s 

4x 10 s r 

2.34 x 10 s 

1.92 

X 

10 s 

2.39 

X 

10 s 

1.96 

X 

10 s 

- 

InN 

9.31X10 1 

1.05xl0 2 

- 

1.06 x 10 2 

1.10 

X 

10 2 

1.14 

X 

10 2 

1.16 

X 

10 2 

- 

InP 

7.38X10 1 

8.52x10 s 

9.3x 10 s [ 

] 6.78 x 10 s 

6.49 

X 

10 s 

6.93 

X 

10 s 

6.82 

X 

10 s 

- 

InAs 

2.04X10 1 

2.52x10 s 

3x10 s a 

1.83 x 10 s 

1.70 

X 

10 s 

1.85 

X 

10 s 

1.80 

X 

10 s 

- 

InSb 

1.26X10 1 

1.38x10 s 

2x10 s a 

9.64 x 10° 

9.56 

X 

10° 

9.76 

X 

10° 

9.86 

X 

10° 

- 


BeO 

3.50X10 2 

4.47xl0 2 

- 

2.80 x 10 2 

2.93 

X 

10 2 

3.27 

X 

10 2 

3.14 

X 

10 2 

3.7x10 2 [ 

] 

BeS 

1.43xl0 2 

1.57xl0 2 

- 

1.19 x 10 2 

1.12 

X 

10 2 

1.26 

X 

10 2 

1.18 

X 

10 2 

- 


BeSe 

9.50X10 1 

l.OSxlO 2 

- 

7.57 x 10 s 

8.02 

X 

10 s 

8.19 

X 

10 s 

8.31 

X 

10 s 

- 

BeTe 

7.75 xlO 1 

9.04x10 s 

- 

6.41 x 10 s 

6.07 

X 

10 s 

7.10 

X 

10 s 

6.29 

X 

10 s 

- 

MgTe 

1.48X10 1 

1.47x10 s 

- 

1.14 x 10 s 

1.27 

X 

10 s 

1.16 

X 

10 s 

1.27 

X 

10 s 

- 


ZnO 

5.04X10 1 

6.45x10 s 

- 

4.10 x 10 s 

5.56 

X 

10 s 

4.13 

X 

10 s 

5.73 

X 

10 s 

6x 10 s 


ZiiS 

3-lOxlO 1 

3.65x10 s 

2.7x10 s 

* 3.15 x 10 s 

3.40 

X 

10 s 

3.24 

X 

10 s 

3.83 

X 

10 s 

- 


ZiiSe 

1.49X10 1 

1.56x10 s 

1.9x10 s 

[ 1.37 x 10 s 

1.42 

X 

10 s 

1.38 

X 

10 s 

1.45 

X 

10 s 

- 

ZnTe 

1.84X10 1 

1.89x10 s 

1.8x10 s 

1.38 x 10 s 

1.48 

X 

10 s 

1.38 

X 

10 s 

1.51 

X 

10 s 

- 


CdS 

1.90X10 1 

1.99x10 s 

- 

1.61 x 10 s 

1.83 

X 

10 s 

1.66 

X 

10 s 

1.85 

X 

10 s 

1.6xlO s [ 

] 

CdSe 

1.03X10 1 

1.14x10 s 

- 

8.42 x 10° 

9.16 

X 

10° 

8.37 

X 

10° 

9.87 

X 

10° 

- 


CdTe 

6.28x10° 

6.67x10° 

7.5xl0°[ 

] 4.90 x 10° 

5.52 

X 

10° 

4.87 

X 

10° 

5.71 

X 

10° 

- 


CuCl 

1.44x10° 

1.26x10° 

- 

9.95 x 10~ s 

1.71 

X 

10° 

1.02 

X 

10° 

1.97 

X 

10° 

8.4x 10~ s 

C 

CuBr 

2.52x10° 

2.75x10° 

- 

1.76 x 10° 

2.83 

X 

10° 

1.74 

X 

10° 

2.99 

X 

10° 

1.25x10° 

c 

Cul 

6.45x10° 

7.10x10° 

- 

5.82 x 10° 

6.85 

X 

10° 

5.99 

X 

10° 

7.22 

X 

10° 

1.68x10° 

c 

CuH 

1.39X10 1 

1.44x10 s 

- 

7.81 x 10° 

6.16 

X 

10° 

8.22 

X 

10° 

7.60 

X 

10° 

- 


Agl 

2.01x10° 

2.44x10° 

- 

1.38 x 10° 

1.78 

X 

10° 

1.37 

X 

10° 

1.89 

X 

10° 

1.03x 10°J 

a 

SiC 

3.55X10 2 

3.72x 10 2 

- 

3.55 x 10 2 

3.11 

X 

10 2 

3.67 

X 

10 2 

3.14 

X 

10 2 

4.9x 10 2 p 



a Ref. [39] 
b Ref. [IQ] 
c Ref. EH 

d Ref. E2I 
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